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ABSTRACT 



The measurement of the human descnihing turret ron and' 
remnant in a compensatory tracking task, la undertaken.. 

These measurements are obtained through the. application, of-' 
the fast Fourier transform technique on a. hybrid (anarog-— 
digital) computer. This method processea the: data:^ in. reall 
time with minimal core storage and the results: are: acvaHable' 
immediately upon completion of the trac]cing run.. 
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I. 



INTRODUCTION 



A. BACKGROUND 

Many of the important tasks performed, by pilots' are? edcin 
to those performed by linear servomechanisms .. Tn. altuaflons: 
such as this, the pilot can be modeled by a set. of: constant — 
coefficient linear dif feren-tial equations.. In -the- frequency 
domain, such a model is often referred -to as a ''h-uman. pilot: 
describing function." The term "describing function" is? 
preferred to "transfer function" -to emphasize the fact 'that 
this pilot model is approximating a nonlinear element, and. 
is valid only for the particular inputs,, system dynamics' and 
task at hand. 

The pilot-describing function is useful in s'tud^inq' two 
classes of problems. First, the describing functions 
measured in the piloted simulation of a. given aircraft. 
and task can be utilized in the subsequent stability and 
control analysis of this aircraft. Once the pilot's 
describing function for a particular task has been measured, 
he can be analytically replaced by his describing function 
in the analyses normally associated with the study of linear 
feedback systems. Second, actual measurement of pilot 
describing functions in ground simulation or in', flight tests? 
can be used to determine how a particular aircraf-t. or f3_ight. 
task affects pilot behavior. Knowledge of pilot describing 
functions consequently provides valuable information for 
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the aircraft designer. 

It is the problem of human describing function. measure- 
ment which forms the basis of this research... The hYbridl 
computer program which has resulted can. be, utilized in. 
virtually any study involving human behavior' in. compensatory 
tracking tasks. 

B. COMPENSATORY TASKS AND QDAS I-LINEARI ZATION 

The compensatory tracking task:,, shown, in: Eigure 1,, 
assumes that the error signal is the only information that 
the operator is receiving. In this study the. operator 
attempts to minimize a visual error signal by using: a hand- 
operated controller. Tracking situations such as this, are 
often encountered in aircraft flight control';; e.g..,, a pilot 
attempting to maintain some desired pitch attitude in. the 
presence of atmospheric ttirbulence .. It has been, shown that, 
in tasks such as this, the operator is nonlinear and' time 
variant in behavior. He may, however, be successfully 
modeled in a quasi-linear fashion [Ref. I] .. This quasi- 
linearization implies that his response ta visual stimula- 
tion is largely linear and time invariant; i..e. , his 
dynamics are largely those of constant-coefficient linear 
differential equations. To account for nonlinear and/or 
time-varying behavior, the model includes a remnant signal 
as shown in Figure 2. The remnant is that portion, of: the 
operator's output which is not linearly correlated: with the 
input. The human operator model thus consists of a linear 

describing function, Y ( jo)) , determined from the quasi-linear 

P 
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analysis, and a remnant, n(t). 

It should be noted that this quasi-Iiaieair modell Lsi ofr 
little use if the remnant is relativ^Xy Iarg_e,, s±ncsi then', 
the operator's behavior is predominan-tLy nonlinear.' and/on.' 
time varying. 

To determine the human operator model i.t is; necessary 
to calculate the linear describing function (j.w) and.' the 
remnant power spectral density, r from- physioall 

measurements of signals of finite duration in. a. laboratory ■ 
experiment. The input signal must appear' to the operator' 
to be random, although it need not be truly random,, andl 
the operator must be well trained;- i..e... ,, not. undergoing 
adaptation or learning [Ref. 21 . 

In order to measure the human describing function. and 
remnant, one of -three techniques can be employed.. These 
are the direct Fourier analysis of the system signals,, the 
use of crosscorrelation me-thods, and a model optimization 
technique. These methods are discussed in Ref.. 3,.. 

If, as done here, direct Fourier analysis or cross- 
correlation methods are employed in measuring the human 
describing function and remnant, then the concep-t of' 
spectral analysis must be introduced. The next section is 
devoted to this topic. A more thorough treatmen-hcan.be 
found in Ref. 4, Appendix D. 
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II. SPECTRAL ANALYSIS 



A. PERIODIC SIGNALS 

A periodic signal, f (t) , witii a fundamental ffequencY 
and period T, satisfying the Dirichlet candi.t±ons' [Ref.. 4.,, 
p. 579] , may be represented by a Fourier series 



where 



f(t) =l F(n)e^^l^ 

n=-oo 

(n) = -^/*'^/^f (t)e‘^'^l"dt.. 

^ J-T/2 



The autocorrelation function for the above periodic 
signal is defined as 



(t) 



= 

y_T/2 



f (t)f (t+x)dt.. 



This can be written 



(t) 



= —f 

Tj. 



T/2 n=-oo 



and is equivalent to 



♦ff(T) = 



1 F(n)e3"“lW4-) 


f ^/^f (t)e^^‘\^dt 


[n=-» ' ^ ’ 


J -T/2 



With F (n) denoting the complex conjugate of F (n) , 

4>ff(T) = I F(n)e^'"“l'^ F(n) = I F(.n) f (n.) 

n 

oo 

'J’ff (T) = i; |F(n) 

n=-oo I 



n=-oo 



n““-“0o 



Zgjnoij^T 
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•i 






The power spectral density, $^^(n) ,, is defined 



$ff(n) 



^ y-T/2 



and it can be shown that 






F(n) 



It can be seen from this relationship that 



n=— <= 

The crosscorrelation function, ^ (x) , of two: periodic 

12 

signals may be found in a similiar manner if both, signals 
have equal fundamental frequencies,, w^,, and both, signals, 
satisfy the Dirichlet conditions^ Assuming- these con dif ions 
are met, then 

ijr y' 2. 

4>f f (X) = f^(t)f2,(t+T)dt .. 



In a fashion similar to that for single signals, it can be' 
shown 

({Jf f (x)= I (n)F 2 (n)e^^^l'^ 

1 2 n=— “ 

The cross-power spectral density, - (n) ,, is defined 

^1^2 



$ 



(n) = r 

12 ^ J-' 



T/2, , » -inco,x,. 

<i>.p f (X) e.-* 1 dx 

T/2 ^1^2 



and it can be shown that 
Using this relationship 



4>f^f^(x) = I ^ (n)e^^“l^ 



f , f- 
n=-“ 1 2 
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B. TRANSIENT SIGNALS 



A signal, f (t> , is defined to be txansient if: 

Lira, f ( t) = Q. .. 



If the transient signal, f(t), satisfies the Dirichlei:. 
conditions in any finite interval, and if 




f (t). 



dt < 



00 



ft 



then the signal may be expressed as a Eourierr' integxaH [Ref. 5', 
p. 279]. Under these conditions, the: Fourier.* integral'. 



f (t> 



2ir / • 



F ( j co) e^'^^dcu 



gives the values of f (t) at all points, including those where the. 
function is not continuous. The Fourier transform of' f(t) is.- 

y^OO 

F(joj) =/ f(t)e ^‘^^dt 



The autocorrelation function for the nonperiodic signal 
is defined 

f 00 

f(t)f(t+T)dt 






/ 

— o 



This can also be written 






= —( 

2tt / 

J — ' 



F ( jcu) 



2_ja)T, 
e-" dto 



Letting 



Cu)) 




/ ! 



where (w) is defined to be the energy spectral density of 
the signal f (t) , it can be shown that 
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■/ 






Thus it can be seen that the energy, density spechrura and: 
the autocorrelation function of a transient signal! are_ aa 
Fourier transform pair. 



♦ fftx) = 






<Ef^(a)) e^^^du) 



and 



-r 



^ -oo 



If two transient signals, f^(t) and f 2 ^(t) each satisfy 
the Dirichlet conditions in all finite intervals, , and. if' 



/ 



fj^(t) 



dt < 0° 



and 



f2<t) 



dt < “ 



then 



/• oo ^co 

4>f £ (t) =/ f (t) f^ (t+T) dt = •^^/ Fj^(lai)F2-(xco) 

12 J 



dcj . 



Now with 



$ 



^ ^ (o)) = F ( jo))F ( jw) , 
12 ^ 



where ^>.p - (oj) is defined as the cross-energy spectral. 

^1^2 

density of the signals f ^ (t) and f^ (t) , it can be shown that 



: f (w) =f (p ^ (T)e^“^dT ,, 

■12 y-oo ^1^2 



and has as its inverse transform 



^fj^f2{T) = IF 






(u) e^^^dai 



The transient signals fj^(t) and f‘ 2 -(t) are said’ tO' be 



linearly uncorrelated when (p^ ^ (x) = 0 for all x. 

• ■ ^1^2 
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C. RANDOM SIGNALS 



In general, a random, signal,, f(t). „ froin a. stationary,, 
ergodic random process [Ref. 4^ p Z73] ,, does not:: have, a 
Fourier transform since 




f (t) 



dt 



is not finite. An autocorrelation function, may be. defined, 
for the random signal f (t) as 




(T) 



LIM 




f ct) f ( t+r) dt .. 



Since (t) satisfies the Dirichlet conditions for: all', 
finite intervals and 




♦ fjCT! 



dr 



< oo 






it can be represented by a Fourier integral.. It: cair: then, be 
shown that 

^ — CO 

where (u) is defined as the power spectral density of' 
the signal f (t) . is the Fourier transform of. 

'fffCw) =1 c|3^^(T)e ^‘^^dx 

If there exist two random signals, f^ (t) and f 2 ^ (t), 

which are sample functions from two different random processes, 

each of which are stationary, ergodic,. and jointly ergodic , 

then the crosscorrelation function is given, as 

4>f f (T) = LIM — (t)f, (ttT.)'dt =- ^/ .. (w)e^'^'^du), . 

12 T->« 2T J-T ^ ^ J-«, ^1^2. 
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where ^ (to) is defined as tiie crass-powex specfraX 
^12 

density of the random signals fj^(t) and f^(t).. ({) ^ ^ ( t) has 

the Fourier transform 



^ (OJ)’ 

12 






The random signals f^(t) and are-, said' to: be: 

linearly uncorrelated if ^ (x) = CT for- alX x:.. 

12 
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III. DESCRIBING FUNCTION AND REMNANT RELATIONS 



A. FREQUENCY DOMAIN RELATIONS 

Again consider Figure 2 witli the ingut e sample:: fiinctioir 
from an ergodic random process. It is seen, that: 

ECjo)) = Kjw) - C(i(jj) , 

where the Fourier transforms are as defined: in Section R-, . NOw 
C(jo)) = 



N(juj) +- Y (ju)E('iu) 

Ir 






then 



E(ju)) = 



Kju) - N(ja))Y^(jo)) 



1 + Yp(jo3)Y^(jo) 



Finally, after multiplying by I(jco) ,. the complex, conjugate 
of the input. 



I(jcj)E(jw) = 



I ( jo) 



I (jo) - N(jo)Y^(jo) 



1 + Yp(jo)Y^(jo) 



In a similar manner. 



P(jo) = Y (jo)E(jo) + N(jo) 

ir 



and then substituting for E(ju) from above, 
[Yp{ju)I(jo) - N ( jo) Yp ( jo) Y^ ( jo) 



P (jo) 



1 + Y ( j o ) Y ( j o ) 
P ^ 



-t N( jo) 



and 



_ Kjo) [y (jo)I(jo) + N(jo) 

Kjo)P(jo) = 



[1 + Y ( jo) Y ( jo) 



p'^ ' e 

Likewise, I(jo)C(jo), E(jo)P(jo), E( jo) E ( jo) ,■ E ( jm) P ( jo) , 
and c (jo)C(jo) may be calculated. The results,, te be 
utilized shortly, are shown in Table I., 
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B. FINITE RUN LENGTH 

Since it is imposstbie to. have experimentaL runs' off 
infinite duration, measurements using finite time histories, 
are necessary. Reference indicates: that: finiter run-, 
lengths can be handled analyti'caliy as- follows:; Iff i^^-^t) 
is the input and defined 

(i i (t) ,, - T £ fc t; 
i^Ct): = , 

^ a „ ELSEWHERE. 

then this function can be considered’ to. be transient: and. 
have a Fourier transform 



iCj'w) = 






^lait 



dt 



The other system signals and their: transforms can be., defined’ 
in precisely the same manner.. It the run: time, T.”,, is. large. 
enough to ensure accurate power spectral, measxirements , yet ^ 
finite so that the respective Fourier- transforms- exist, 
then the following spectral relations are valid [Ref... 6] . 



and 



* / ^ 



♦ ■ M = ™ 



2T 



I( jm) 



2^- 1( jo))P ( joj) 



where and power; and cross-power.- spectral 



densities respectively. Utilizing Table I, 
LIM \ _1 



$ . (oj) 
xp ' ' 



but 






I 



2T 



I(jw)I(jia))y. (joj) -t I(jtjj)N(ja))Y (jw) 
P ^ 



1 + Y (jto)Y (jto) 



LIM 



2T 



I(jo))'I(i(jj) 



— $ . . (cj) 
rx 
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and 



Lm 



I( jai).N.(i(i3) 



= (OJ) 

xn 



thus 



Y C joi): <5 . .. (o) ) +- $ . ( 0) ) Y_ { jto ) . 

/Fv / \ _ T !■ iXl O 

^ip “ 1 -h Y ( jcj) Y ( JO)) 



Now, 



= jT <^^.^{(jj)e 



xn 



xn 



Since by definition the rremnant,. n.(t) ,, ils- rineaxry.' uncorrelated _ 



with the input, i (t) , then 2 for all t:.. Thuss 



$. (u) = 0 and 

xn 



Y (lo)) $ . . (u) 
-K / \ p XX 

$. {(ji) = 

ip . 



1+ Yp( .ja)) Y^(jto) 



In a like manner. 



* /■ t LIM 

$. Cw) = „ 
xe T-»-<=° 



L(iu)E(jco) 



or 



(w) = 
xe 



. (to) - (to)Y (ito) 

XX xn c -* 

1 + Y (j(j)Y (j(j) 

p ' j ^ ' j. 



Again since n(t) and i(t) are linearLy uncorrelated. 



= 0 for all T, then =0.. Thus~ 



$ •. (to) — 
xe ' 



. (to) 

XX 



L+Yp(jto)Y^,(jto) 



Also 



<f. / ^ 'hIM 

$ (to) = _ 
pp T->“ 



P(j(o)F(jto) 



$ (to) 
PP 



LIM 


) 1 


■ ^(jto)Y (j(o)+f7:(jto) 


'r( jTO ) Y_ ( j to) -t-N ( j to ) 


T-><» 


(2T 


-1 + ■ 


-1 + Y (jto)Y (jto) - 



I' 
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or 



$ ((jj) 

PP 



^..(ti>)iY { joj) I (joj)4-$. (oj)Y (joj)+4 (oj) 

11 1 p 1 ai p in P nn' 



L+r (jw)Y^.(jw) 
p c. • 



Again, since 



« . (w) = Q,. 
ni 



$ {(ij) 

PP 



4-- (w) 


Y Cjo)) 


+4 (to) 


11 


-P_ 


nn 



L+Y (joj)Y' (joj) 

p\j (,\j 



C. DESCRIBING FUNCTION AND REMNANT 



From (co), (co) , and $ (w) ,• Y' ( joj) and: <1 (w) may 

ip xe pp. p nn. 

be found. Utilizing (w) and $. (to). 



iP 



$ii((o)- = [l+Yp( joj) Y^ (i(o)| 



and 



11 Y (ito) ^ ^ 

P ■' 



Thus, 






or 



«.„(C0) 

Y (joj) = ^ 



Also, from $ (to) , 
PP 






$ (to)- 

. PP 



Yp(ito) 






In addition, Y^Cjjo) can be determined and calculated from 






Again = 0, thus 
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^ • { (u) 
xc 



xe 



(ui) ( joj),Y^ ( j;w) 



But 



('to),T (tw) = $:. (w) 

xe p xg 



thus 






■t. 



XC 




(to) 



The fxxnctions in Table IX form ther banxs: of. the3 describing 
function measurement technique: utiTized' in this* studjr.- . 



D. SINUSOIDAL INPUTS 

The relations in Table IX are predica.ted' on. the; existence ‘ 
of a random input. In experimental work, a* random appearing ^ 
input is often used and can be generated' as a summation of ' 
sine waves [Ref. 31- The input can. thus be represented byy 

n 

i('t); A, sinco. t 

k=l 

where the are chosen such that in. a finite: run-, length 
there will exist an integral number of periods' orr complete 
cycles. In this analysis a run time of 150 seconds was 
used and 0.08 < ^ 40-0 radians per second. 

Utilization of a sinusoidal input results in. system 
signals that have both random and periodic components. With 
what is now a mixed signal, the question, arises as to which 
power spectral relationship should be used.. The, solution, is 
to use the periodic power spectral relationships^ forr measurer - 
ments made at the input frequencies and to use the, random, 
finite relationships at all other frequencies. 
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It should be noted that if the: eiiperimentail mn‘ length, 
T, is large and contains an. integraf number.' off periodss of i 
each of the input sinusoids,, then, the Eourfer. transforms- of f 
the periodic and finite random signais differ* onrly/ b.y/ai 
constant of proportionality [Eef.. 6T- 

In the periodic case, it was shown, that. 

F^(n) = „ 

and for the random signal , 



/•oo 

7 



Ffjm) = / ^(t)e 



-jolt 



dt 



If it is recalled that f^(t) = 0. for- t<-T and t>T:,, and_ 
letting T = m^^Tj^, where m^^ is the niiraber' of periods, Tj^, of 
frequency then 



k^-T^ kk^^-iu^T^^ 



or 

AT 

Fj^(n) = ^f (t)e"^“k^dt ;; 

thus it can be seen that 

2TFj^(n) = E(jo)j^) .. 

It should be noted that the same expansion', off the; limits 
on the integral can be used for the periodic, function; thus , 



$ 



ff 



'"k> 





f (t) f (t+T) dte 



-3^“k^d.T 



t 
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or 






' a/V’ 

4T ^ -T J -n 



f (t) ECt+T) e 



4T -T -T 
Also for the finite random function,, 

'T 






$ 



ff 



(w)^ = f^.(t) f^(ttT.)e~^’^dtdT: ;; 



or 

►T /rT 

J i 

-T *^-T 

then by equating the integrals,, since f^(t) = f:(t) fbrr 
-T<t<T, it is seen that at a specific- input frequency , oj,, , 






where 



^ _ lim: 

^ff^ k^T 



2T 



F ( ioij^ ) 



It has been shown that 



$. (w) 

Y (jw) = 



(cu) 

xe 



At the input frequencies , with T large and containing an 
integral number of periods of each input frequency, 

<5. (o), ) <l> . (w, ) $ . (n, ) 2T 

Y fTM ■» = J-P k ^ iP k T ^ iP ^ 

k^ «. (o, ) 4>. (oj, (n, )2T 

^ xe k xe' k'T xe k: 



or 



Y (joj ) 
P k 



*ie'V 



L(n,J P (n-j^) 



F(n^) 






Similarly , 
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*ic'V C'V 

♦. (n.) P(n^) 

ip ' k k' 

This illustrates that the cross-power spectral measure- 
ments need not be made and only the Fourier transforms are 
needed. The latter is usually an easier measurement than the 
former . 

The terms in $ (w) can be examined in a similar manner, 

nn 

In this case, if measurements are taken at frequencies, 

other than those used in the input, then the expression for 

<l> (to) is somewhat simplified; i.e., 
nn 



$ (w. ) = 

nn h 



1+Y (j(0,)Y (jw, )l^ $ ((0.) , 

p h c h' I pp h ' 



since $ . . (to, ) = 0 . 

11 h 

It should be emphasized that at frequencies, other 

than those used in the input. 



P(nj^) 



C(nh> 

^ P(V 



since I(n, ) = I(n, ) = 0. Thus Y (jto, ) and Y (jto. ) must be 
n n p ii c n 

estimated, as direct calculations can be made only at the 
input frequencies. In order to estimate ^^(jtOj^), simple 
linear interpolation between Yp(jtOj^) and 
used. 
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IV. COMPUTER MECHANIZATION 



A, EXPERIMENTAL SET-UP 

The measurement of a human's describing function and 
remnant in the compensatory tracking task of Figure 2 was 
made using a hybrid (analog-digital) computer. The error 
signal, e(t), was viewed as the vertical displacement of a 
horizontal line on an oscilloscope screen. ‘The operator's 
controller consisted of a non-moving force stick. Control 
was effected by fore and aft pressure on the stick; e.g., 
if the line on the oscilloscope moved above the datum, the 
operator applied forward pressure to move the line down, and 
vice-versa. The input, i(t), and controlled element dynamics, 
Y^(jw), were mechanized on the computer as were the measure- 
ment algorithms to be described. Each experimental tracking 
run lasted 150 seconds. 

B. FAST FOURIER TRANSFORM 

The pertinent relationships are again 




and 





Now 




-T/2 
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s?. 






or 



P (Hj^) ^ 



^T /2 

I p ( t) cos(rui)^t) dt +- jf p:(t) s±n.(nw, t) dt 

''-T/7 ^ J'-rr/T: 



T/Z 



T/Z 



L*'-T /2 

These integrals may be approximated, by the. follow ingr 
summations : 

P(n.) = 2 p(nAt)cos ((Oj^nA_t)+i ^ ]^'p (nAt) sin (u^nAt) 

n=0 n=0- 



If 



K 

E 

n.=0 



coe(Wj^.nAt) 



and 



B 



K 



Pj. sin.(a)^nAt) 

n=0 



then the fast Fourier transform E (n, ) can. be. written', as. 

k. 



p ^ ^ 



A + JB 

: Pk. Pk 



Similarly, it can be shown that 
At 



C(n^) = ^ 



A + j B 
c, c, 

k k 



and E.(nj^) =5^; ^ 



.. At A t- j3. 



I 



From this , 



and 



A. +- f B 

^p ^ ^ A +■ j B 

^k- ^k 






It can also be shown that 



-iL 2 B 2 

l?k Pk 

2 2 
A + B 
e, e, 

k' k: 



1/2 



^ = tan' 



- 1 . . -1 Pk 

— taui'. -5 



A. 



®k 



^-k: 
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In order to validate the siimilatedl cantroiled ehenrent- 
dynamics, Y (jw) , on-line measurement o£ Y.(j_ai) can. be:; 

C C- 



utilized; 



and 






L , B; 



2 . 






^ Y^tjojj^) = tan. 



Z 2 

A + E 



r/z 



Bk : 




E 


& 


Bk -11 

A. 


0 




Bk: 


^k: 



The power spectra of the remnant can also: be' determined 
from the above relationships with the interpoiation . process 
described earlier. The determinatton. of the power..' spectra 
of the operator output, accomprished'l using 

the measurements p(t) at any desired frequency [Refi. 6] .. 
Previously it was shown that 

pp h^ pp h T pp h ' 



thus 



• VlJU, J 

PP ^ 



2 



h' ' h^ 





2 : = (^)^2T. 


A ^ -t b: ^ 


h 


T^ 


Ph. Ph. 



and finally, 



. , , ‘ 2 (At) 

(w. ) = — t — '' . 
pp h T 



^ B ' 

B, Pu 

h. h: 



From this it can be seen that in order to determine the 



operator's describing function, Y' (301,),. the remnant, 

p K. 

^k! 



/ and the controlled element dynamics,, Y. _(jWj.) , the 



only measurements needed are those of the' error, opemtor. 
output, and controlled element output. rf each of: those 
measurements, taken at specific times, is then multiplied 
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by the proper trigonometric function, and summed over' the: 
entire run , then the describing functions: and remnant can . 
be calculated. 

The major drawback of the fast Fourier' transform 
technique is the necessity of calling the trigonometric:: 
functions sin(x) and cos (x) during the rtm... This requires 
so much computer time that the multiplication, and" addition 
computations cannot be performed between. analog'-te-digitaU 
interrupts. This in turn means that. alL data, must be stored 
for later computational purposes. 

The necessity of calling trigonometric functions during 
the run can be avoided if the relationships for' sin (a+b) and, 
cos (a+b) are utilized in recursive fashion.. If the initial 
time of the run is "a" and the time between measurements is 
"b" then 

sin (a+b) = sin (a) cos (b) -t- cos (a) sin (b) 

and 

cos (a+b) = cos (a) cos (b) +- sin (a) sin (b) .. 

It is obvious that a recursive process can be mechanized 
which obviates calling trigonometric functions during the 
run. 

The use of this method yields results immediately upon 
completion of the run and conserves computer, storage space. 
Within seconds of run completion,, the. data may be .analyzed 
from either numerical or graphicaX read-out.. 
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V. RESULTS 



The operation of the program was checked by measuring 
the transfer function of known elements or fiLters: iir. 
place of the operator . The results are shown, in . Figure. 3', , 

A, 5, and 6. 

In validating the remnant measurement: technique:,, the. 
operator was again replaced by an. element- with known- transfer, 
function. A signal 

n{t} = Asin (Wj^t) " 

was inserted at the output of this element.. Measrured' 
values of were then compared with, the theoretical; 

j value [Ref. 71. 

Actual describing function and remnant measurements, 
were then taken on two subjects .. The first: oi these, had' 
considerable previous experience.. The second' subject had 
experience only with this experimental set-up.. The measure- 
ment results are shown in Figures 7-18 for controlTed 

elements of Y (s) = l.Q and — . Each of ' these figures 
c s 

represents the results of 10 tracking runs.. In. all runs, 
the human describing function contains the gain of; the 
controller . 

The describing functions illustrated are comparable 
with those obtained by other experimenters;; e..g...,, [Ref:.. 8] .. 

This computer program represents a powerful tool, for 
use in experimental investigations involving a human 
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controller. It can, for example, be utilized in, a variety^ of.' 
situations in v/hich quantitative models of pilot behavior: 
are desired. The program itself req\iires Irtfle- cores storage:. . 
This means that considerable storage is available- for.* simu- 
lating complex aircraft dynamics , providing detailed disp.lay_ 
formats and calculating performance measures.. 
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TABLE 1 



E(jw) 






C(ju^) 

T(jaj)E(ju;) 






T(j u^) C(j o^) 



E(ju>) E(jco) 



P(jaj)P(joj) 



C(ju>)C(jo;) 



Uj‘^) - 

I + Yp(joj)Yc(jcj) 



• (jtJ) Yp(j oi) + N(jcj) 
I + Yp(juJ)Y^(jaj) 



= N(joj)Yj.(j o>) + E(juJ)Vp(ja)) Y^.(jw) 



T(jaj)Kjcj) - T( j u^) N ( j w) Y^(jtJ ) 

i + Yp(joj) Y^(joj) 



Kjoj ) I (jcj ) YpC j ) + T(ja^)N(ja^) 

~ Yp(jc*^ ) Yj.(j cj) 



— [T(joj) N( joi) + T(joj ) E (j oj ) Yp(joj )] Yc(jcu) 



Kju) ) - N( ju ) Yj,( jcu ) 




l(j oj) - N(juj) Y^U <^) 


_» + Yp(j(j) Yp(j u^) 




_ 1 + Yp(j 0J)Yc(jaJ ) 



l(j‘^) Yp(joj) 4 N(j CJ) 


|(ja;)Yp(jaJ ) + N(jcJ) 


_ 1 + Yp{ju»)Yp(j a;) 


1 + Yp(ju;) Yj.(j oj) 



= [N(joj) + E(j w)Yp(jujj|Yj.(jcj)][N(jw) + £(j cj) Yp(ju) ) [y^O oj )] 
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TABLE IC 



Yp(jw) 
Yp(jo;) 









" <l>il(w)IYp(jaj)l^ 

<&ii(^)IYp{jo))l^ + 

11+ Yp(juj)YJjco)!^ 
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HUMAN OPERATOR in a COMPENSATORY TRACKING TASK 
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FIGURE 



QUASI-LINEAR MODEL of HUMAN OPERATOR in a COMPENSATORY TRACKING TASK 
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FIGURE 8 
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AMPLITUDE of HUMAN DESCRIBING FUNCTION for Y^(s) = 
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